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We determine the (pseudo)critical lines of QCD with two degenerate staggered fermions at nonzero 
temperature and quark or isospin density, in the region of imaginary chemical potentials; analytic 
continuation is then used to prolongate to the region of real chemical potentials. We obtain an 
accurate determination of the curvatures at zero chemical potential, quantifying the deviation be- 
tween the case of finite quark and of finite isospin chemical potential. Deviations from a quadratic 
dependence of the pseudocritical lines on the chemical potential are clearly seen in both cases: we 
try different extrapolations and, for the case of nonzero isospin chemical potential, confront them 
with the results of direct Monte Carlo simulations. Finally we find that, as for the finite quark 
density case, an imaginary isospin chemical potential can strengthen the transition till turning it 
into strong first order. 

PACS numbers: 11. 15. Ha, 12.38.Gc, 12.38.Aw 



I. INTRODUCTION 

The determination of the QCD phase diagram in the 
temperature - quark density plane is becoming increas- 
ingly important, due to its impact in cosmology and in 
the physics of compact stars and of heavy-ion collisions. 

The first-principle nonperturbative approach of dis- 
cretizing QCD on a space-time lattice and performing 
numerical Monte Carlo simulations is plagued, at nonzero 
quark chemical potential, by the well-known sign prob- 
lem: the fermionic determinant is complex and the Monte 
Carlo sampling becomes unfeasible. An exact solution to 
the problem is yet not known, but various approximate 
alternatives have been explored. One possibility is to 
investigate other models which are free of the sign prob- 
lem, like two-color QCD or QCD with a finite isospin 
density [l|, [H; however in such approach predictivity is 
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plagued by the fact that systematic differences between 
the model and the original theory are not known a priori. 

A positive measure is obtained also when the quark 
chemical potential is purely imaginary, in this case the 
idea is to infer the behavior at real chemical potential 
by analytic continuation. Such approach was first sug- 
gested in Ref. [3j], while the effectiveness of the method 
of analytic continuation was pushed forward in Ref. 
Since then, the method has been extensively applied to 
QCD with staggered 0, i, 0, i, i, [13, p, [ll and Wilson 
quarks [H, [ij] and tested in QCD-like theories free of 
the sign problem [H [H, [13, [ll [H, [13, [H and in spin 
models [22, 23l. 

The idea underlying the method of analytic continua- 
tion is very simple: if the dependence of an observable or 
of the critical line itself on the imaginary quark chemical 
potential is expressed in terms of an analytic function in- 
side a certain domain, then this analytic function can be 
prolongated to the largest possible domain, compatible 
with the presence of singularities, up to the physically 
relevant region of real chemical potentials. 

There are, however, two important limitations to the 
effectiveness of the method: a practical one, due to the 
fact that Monte Carlo simulations yield data points (with 
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statistical uncertainties) at fixed values of the imaginary 
chemical potential and, therefore, analytic continuation 
passes through the choice of an interpolating function, 
which may be ambiguous; a principle one, due to the non- 
analyticities and the periodicity of the theory with imag- 
inary chemical potential which makes so that the 
region effectively available for Monte Carlo simulations 
is limited by the condition Im(/i)/T < 1. The combina- 
tion of these two drawbacks implies that the region of real 
chemical potentials where the analytic continuation is ex- 
pected to be reliable can be estimated as Re{fi)/T < 1. 

Apart from analytic continuation, a careful study of 
the phase diagram in the T - Im(/i) plane is important 
by its own. It may teach us something about the critical 
properties of QCD also at zero or small real n [1, [H, 
m, 113, m, and, at the same time, it can be used to 
the test the reliability of QCD-like models which are then 
used to explore real fj, as well [13, HH, IH, [H, [ll, [H, HI, 

113, m,!!!. 

An important role is played in this respect by the pe- 
riodic series of unphysical first order lines, located at 
Im(^)/T = {2k + l)TT/Nc, fc = 0,1,2,..., being the 
number of colors, which characterize the high-T region 
of the T - Im(^) plane (Roberge- Weiss lines). Such lines 
are connected to the analytic continuation of the phys- 
ical pseudocritical line by an endpoint which, both for 
rif = 2 and n/ = 3 QCD, is first order (triple point) 
in the limit of small or high quark masses and second 
order for intermediate mass values; it has been conjec- 
tured that this may be directly related to the nature of 
the pha se transition at zero or small real chemical poten- 
tial [H, 113, H, H . 

In a series of studies by some of us, we have started 
a detailed investigation aimed at checking the reliability 
of analytic continuation of the pseudocritical line and at 
extending its range of applicability, by looking for pos- 
sible deviations from the simple linear behavior in /i^, 
Tciii'^) = Tc{0)+Afi'^, which fitted well with earlier stud- 
ies d, Q. To this purpose, the range of Im(/i) values 
included in numerical simulations was extended with re- 
spect to earlier studies, up to reaching the border of the 
first Roberge- Weiss (RW) sector, statistics was consider- 
ably increased and different kind of interpolating func- 
tions were considered. In order to validate the differ- 
ent interpolation options, in some cases QCD-like the- 
ories were adopted (such as SU(2) or SU(3) with finite 
isospin density) which, being free of the sign problem, al- 
lowed the comparison of extrapolations with the results 
of Monte Carlo simulations performed directly at real /i. 
A detailed summary of such investigation is reported in 
Section |lTl the main result emerging from it is that non- 
linear corrections are not negligible, but an unambiguous 
extrapolation to real ^ fails for fi/T ^ 0(1). 

In the present study we consider rif — 2 QCD in pres- 



ence of a quark (^q) or an isospin ^ (/Xiso) chemical poten- 
tial, whose partition function, in the standard staggered 
discretization for fermion fields, reads 

Z,/i,o(r,/i) = J VUe-^'^idctM[fi])i{detM[±^i])i ,{1) 

where the plus/minus sign refers to the quark/isospin 
chemical potential case. So is the lattice gauge action 
and M is the fermion matrix in the standard staggered 
formulation: 

1 ^ 

+ \ n^A (e"'' - e-"^^ ■ (2) 

Here i and j refer to lattice sites, j) is a unit vector on 
the lattice, m.v are the staggered phases, a is the lattice 
spacing and m is the bare quark mass. We shall consider 
a bare quark mass am = 0.05, corresponding to a pion 
mass TO^ ~ 400 MeV. 

The partition function in Eq. ([1} is expressed as 
a functional integral with a positive measure, hence 
suitable for Monte Carlo evaluation, when [i is purely 
imaginary, but also when [l is real for Z-^^^^ alone. Our 
plan is to perform an extensive investigation about the 
location and the nature of the deconfinement transition 
for all cases in which Monte Carlo simulations are 
available. The specific purposes that we have in mind 
are the following: 

1) verify the reliability of analytic continuation of the 
critical line, Tc(/i), from imaginary to real /i in the case 
of a finite /liso, where simulations are available both for 
imaginary and real ^iso- Apply analytic continuation to 
the case of a finite /iq, also on the basis of what learned 
in the case of a finite \i.\sa\ 

2) make a careful comparison between the two theories 
at finite /ig or /iiso, quantifying systematic differences for 
quantities like the curvature of the pseudocritical line at 
zero chemical potential; 

3) determine how the nature of the transition changes 
as a function of the chemical potentials. In particular, 
in the case of a finite /i^, no critical point is expected 
on the imaginary side since the adopted quark mass, 
am — 0.05, is slightly above the lower tricritical quark 
mass determined for = 2 in Ref. [l^, hence the 
endpoint of the RW line is second order. The situation 
can be different, hence potentially more interesting, in 
the case of a finite /ijso, where the available range of 



Notice that, in order to permit a direct comparison with /^q, 
we define /iiso = (Mu — Md)/2, where /Xu and [id are the u and d 
quark chemical potentials, i.e. a factor 2 lower than usual. 
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imaginary values is larger. 

Most numerical simulations have been performed on 
a 16^ X 4 lattice; different spatial sizes have been taken 
into account to investigate the critical behavior in a few 
specific cases. The paper is organized as follows. In Sec- 
tion In] we discuss our determination of the critical line 
and its analytic continuation to real chemical potentials. 
Section IIIII is dedicated to a systematic comparison be- 
tween the curvatures of the critical lines at finite fj,q and 
/iiso respectively. In Section ITVl we discuss how the order 
of the transition changes as a function of the chemical 
potential. In Section |V] we summarize our results and 
draw our conclusions. A partial account of our findings 
has been reported in Ref. (40l |. 



II. ANALYTIC CONTINUATION OF THE 
PSEUDOCRITICAL LINE 

Before presenting results for n/ = 2 QCD, it is worth 
making a short summary of our previous findings. 

1) In SU(2) with Uf = 8 staggered fermions and 
finite quark density it was found that the analytic 
continuation of physical observables is improved if ratios 
of polynomials (or Pade approximants [4ll |) are used as 
interpolating functions J 7.] . 

2) In SU(2) with rif = 8 staggered fermions and 
finite quark density [la, l2l[ and in SU(3) with Uf = 8 
staggered fermions and finite isospin density [2l| it was 
found that the nonlinear terms in the dependence of the 
pseudocritical coupling /3c on /i~ in general cannot be 
neglected and that the extrapolation to real ^ may be 
wrong otherwise. Moreover, the coefficients of a Taylor 
expansion in /i^ of l3c{^J?) were found to be all negative, 
implying subtle cancellations of nonlinear terms at imag- 
inary PL in the first RW sector, hence a practical difficulty 
in the detection of such terms from simulations at /i^ < 
only. It was realized that, in general, a 3-parameter 
fit (e.g. an even polynomial of order /x®, with the 
coefficient of the fj? term possibly constrained by a fit 
restricted to smaller /j,^ < values) provided a very good 
description of the pseudocritical line in all explored cases. 

3) In SU(3) with n/ = 4 staggered fermions and 
finite quark density [12 deviations in the pseudocritical 
line from the linear behavior in /i^ for larger absolute 
values of /x^ were clearly seen. However, it turned out 
that several kinds of functions were able to interpolate 
them, leading to extrapolations to real /i which start 
disagreeing from each other for /j,/T > 0.6. In this 
case, contrary to the studies mentioned above, direct 
simulations at real chemical potentials are not available: 
one is not able to decide which extrapolation is the right 
one and the disagreement is in fact a measure of the 
systematic ambiguity related to analytic continuation. 



In the present study we approach the case of SU(3) 
with Uf = 2 and a standard staggered fermion discretiza- 
tion. In principle one expects that issues related to an- 
alytic continuation of the critical line may depend on 
the number of flavors, since the coefficients of the Taylor 
expansion in /i^ themselves have such dependence: for 
instance it is known that the curvature of the pseudo- 
critical line at /i = is smaller for smaller ny, hence the 
sensitivity to nonlinear terms in /i^ could be enhanced. 

Moreover we shall take into account also a finite ^iso, 
considering both the imaginary and the real potential 
case. That will give us the opportunity to directly check 
the validity of analytic continuation and to have a test- 
ground available for the different extrapolations in a the- 
ory which is free of the sign problem and is as close as 
possible to the one explored at finite pq. In view of that 
we shall discuss results at finite fiiso at first. 

The range of imaginary chemical potentials which are 
useful to analytic continuation is limited by the period- 
icity in Im(/i)/T, which is 2tt/Nc for fiq and 27r for ^iso, 
(see for instance Ref. [ll| for a discussion on this point) 
and by the presence of unphysical phase transitions in 
the high-T region. In the explored Nc = 3 case, nu- 
merical simulations will be limited, in the finite quark 
chemical potential case, to Im(/Ltg)/T < 7r/3, where the 
first RW transition line is met at which, in the high-T 
region, the Polyakov loop switches from one sector to 
the other. At finite isospin chemical potential instead we 
limit simulations to lm{fiiso)/T < 7r/2, where at high-T 
a RW like transition is met at which G-parity is sponta- 
neously broken and the Polyakov loop develops an imag- 
inary part [2l|. 

We have adopted a Rational Hybrid Monte Carlo 
(RHMC) algorithm, properly modified for the inclusion 
of quark/isospin chemical potential: modifications with 
respect to /i = arc trivial apart from the case of a real 
isospin chemical potential, where the usual even-odd fac- 
torization trick does not work and an additional square 
root of the determinant is needed. Typical statistics have 
been around 10k trajectories of 1 Molecular Dynamics 
unit for each run, growing up to 100k trajectories for 4- 
5 /3 values around the pseudocritical point, for each /j,^, 
in order to correctly sample the critical behavior at the 
transition. 

The pseudocritical /?(^^) has been determined as the 
value for which the susceptibility of the (real part of the) 
Polyakov loop exhibits a peak. To precisely localize the 
peak, a Lorentzian interpolation has been used. We have 
verified that the determinations are consistent if the sus- 
ceptibility of a different observable, such as the quark 
condensate or the quark number (isospin charge) is used. 
In subsections III Al and III Bl ^ will stand respectively for 

Miso 

and ^iq. 
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A. Nonzero isospin chemical potential 

In Table H] and in Fig. [T] we present our determina- 
tions of the pseudocritical couplings, both for negative 
and positive /i^. As a preliminary step, we have tried if 
an analytic function of exists, able to reproduce all 
the available data, both at negative and positive /i^. It 
turned out that no even polynomials in /i^ up to the 4th 
order can do the job and that the first successful global fit 
is achieved with a ratio of a 4th to 2nd order polynomial 
(see Fig. [1] and Table |TT] for the fit parameters and their 
uncertainties). Also the fit with a 6th-order polynomial 
and the "physical" fit defined below give a reasonable 
global fit, with xVd.o.f. < 2. 

However, it is interesting to notice that a simple lin- 
ear function in /i^ fits well if one includes all data but 
those with < -0.375^ (see 2nd row of Table[IT|. That 
means that, contrary to what we observed in our previ- 
ous studies, in this case non-linear corrections are more 
important for imaginary values of ^ than for real ones, 
where instead, in the range explored in the present study 
and within errors, they are negligible. 



TABLE I: Summary of the values of j3c{y?) for finite isospin 
SU(3) with n/ = 2 on the 16'^ x 4 lattice with fermionic mass 
am=0.05. 



/i/(^T) 




0.475i 


5.41670(31) 


0.4625i 


5.40948(40) 


0.450i 


5.40429(51) 


0.435i 


5.39780(59) 


0.4175i 


5.39012(49) 


0.400i 


5.38353(44) 


0.375i 


5.37588(61) 


0.350i 


5.36799(62) 


0.327i 


5.36239(64) 


0.300i 


5.35570(50) 


0.260i 


5.34820(47) 


0.230i 


5.3425(10) 


0.200i 


5.33800(52) 


0.165i 


5.33304(85) 


0.120i 


5.3289(12) 


0. 


5.32371(86) 


0.050 


5.3199(24) 


0.100 


5.3189(22) 


0.150 


5.31486(82) 


0.200 


5.3091(13) 


0.250 


5.3022(28) 


0.300 


5.2928(15) 


0.350 


5.2788(15) 


0.400 


5.2657(18) 


0.425 


5.26079(94) 



We have tried several kind of interpolations of the pseu- 
docritical couplings at < 0. At first, we have consid- 
ered interpolations with polynomials up to order //^ (see 
Table HIl 3rd to 6th rows, for a summary of the resulting 
fit parameters). We can see that data at /Lt^ < are pre- 
cise enough to be sensitive to terms beyond the order /i^; 
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• lattice 16^ X 4 
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FIG. 1: Pseudocritical couphngs obtained in finite isospin 
SU(3) with n/ = 2 on a 16^ X 4 lattice with am=0.05, both 
at /i^ < (diamonds) and /i^ > (circles). The dashed line 
represents a global fit to all data with the ratio of a 4th- to 
2nd-order polynomial; the solid lines around the fitting curve 
delimit the 95% confidence level (CL) band. 



indeed, a good /d.o.i. is not achieved before including 
terms up to the order /x^. 

As in Ref. 2l[, we performed a "constrained" fit: first, 
the largest interval [(m/(7I'T))^jj^, 0] was identified where 
data could be interpolated by a first order polynomial 
in (/i/(7rT))^, with a y^/d.o.i ^ 1; it turned out that 



(/i/(7rT))^i„ = -0.375^. Then, all available data were 
fitted by a 6th-order polynomial, with the constant term 
and the quadratic coefficient fixed at 5.3232 and —0.368, 
respectively (see Table HIl 7th row). 

Then, we have considered interpolations with ratios of 
polynomials of order up to (jj. / (ttT))^ . The interpolation 
with the least number of parameters for which we got a 
good fit is the ratio of a 4th- to 2nd-order polynomial, 
see Table HJ 8th row and Fig. [^left). 

Finally, we have tried here the fit strategy first sug- 
gested in Ref. [2l| , consisting in writing the interpolating 
function in physical units and to deduce from it the func- 
tional dependence of f3c on /i^, after establishing a suit- 
able correspondence between physical and lattice units. 
The natural, dimensionless variables of our theory are 
T/Tc{0), where Tc(0) is the pseudocritical temperature at 
zero chemical potential, and /j,/(7rT). The ratio T/Tc{0) 
is deduced from the relation T = 1 / {Nta{/3)) , where Nt is 
the number of lattice sites in the temporal direction and 
a(/3) is the lattice spacing at a given /3. Strictly speaking 
the lattice spacing depends also on the bare quark mass, 
however in the following evaluation, which is only based 
on the perturbative 2-loop expression of a{/3) for Nc ~ 3 
and Uf = 2, we shall neglect such dependence. 

We considered the following "physical" fit ratio {x = 
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TABLE II: Parameters of the fits to the pseudocritical couplings in finite isospin SU(3) with Uf — 2 on a, 16^ x 4 lattice 
with fermionic mass am= 0.05, according to the fit function ^c(m^) = (ao + ai(^i/(7rr))^ + a2{tJ./{-!TT))'^ + ai[^ / (TiT)f) / {1 + 
a4{fi/(TTT))^). Blanlc columns stand for terms not included in the fit. The asterisk denotes a constrained parameter. Fits are 
performed in the interval [M/('''"2^))miii! M/(''"r))max]) the last two columns give the value of {l-i-/{''rT))^-^^^^^^. 

ao ai a-2 as ai x^/d-o.f. W(^r))min (M/(7rr))max 



5.32326(62) 16.755(10) -1.072(26) 3.2143(19) 060 -0.475^ 0.425^ 

5.32385(54) -0.3597(60) 0.96 -0.375^ 0.425^ 

5.31940(76) -0.4192(47) 18.3 -0.475^ 

5.3232(11) -0.368(12) 0.59 -0.375^ 

5.3255(14) -0.286(25) 0.511(94) 1.85 -0.475^ 

5.3235(21) -0.374(68) -0.36(63) -2.4(1.7) 0.43 -0.475^ 

5.3232* -0.368* -0.253(91) -2.01(44) 0.62 -0.475^ 

5.32403(94) 14.602(14) -0.844(44) 2.8066(25) 0.49 -0.475^ 




FIG. 2: Fits to the pseudocritical couplings in finite isospin SU(3) with n/ = 2 on a 16^ x 4 lattice with fermionic mass 
am=0.05: ratio of a 4th- to 2nd-order polynomial (left) and "physical" fit according to the function (|4]) (right). 
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leading to the following implicit relation between /3c and 

a2(/3c(M'))|2-ioop = «'(/3c(0))|2-ioopii^:^^^ , (4) 

with these resulting parameters 

^c(O) = 5.32422(94) , A = 4.077(23) , 

B = 2.659(77) , C = 3.221(26) , (5) 

and x^/d.o.f.=0.53. In Fig. [21^ right) we compare the 
"physical" fit ratio to data for /3c (/x^). 

In Fig. [3] we have plotted the extrapolations to the 
interval < ijl/{tiT) < 0.5 of the following fits: 

• 6th-order constrained polynomial (7th row in Ta- 
ble^; 



• ratio (4,2) of polynomials (last row in Table lll| : 

• "physical" fit ratio, Eqs. 

The three curves agree as long as /^/(ttT) < 0.2, but 
then the 6th-order constrained polynomial deviates from 
the other two curves. This means that different inter- 
polations, which all reproduce the trend of data in the 
fit region -0.475^ < (/^/(TrT))^ < and take correctly 
into account the deviation from the quadratic behavior in 
that region, lead to distinct extrapolations, as it occurred 
in 71/ =4 SU(3). We can see that ratio of polynomials 
(Fade approximants) in general tend to be closer to direct 
determinations of the pseudocritical couplings, which are 
reported in the same figure for a few values of (/i/ (ttT))^. 



B. Nonzero quark chemical potential 

In Table IIIII we summarize our determinations of the 
pseudocritical couplings. We have tried several kinds of 
interpolation of the pseudocritical couplings at < 0. 
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FIG. 3: Extrapolation to real isospin chemical potentials 
of the 6th-order constrained, ratio (4,2) of polynomials and 
"physical" ratio fits (only the borders of the 95% CL band 
have been reported). Data points (circles) are the results of 
Monte Carlo simulations performed directly at real isospin 
chemical potential. 



TABLE III: Summary of the values of jScip?) for finite density 
SU(3) with 71/ = 2 on the 16'' x 4 lattice with fermionic mass 
am=0.05. 



Im(M)/(^T) 




0. 


5.32371(86) 


0.100 


5.3277(12) 


0.180 


5.33524(71) 


0.200 


5.33914(83) 


0.245 


5.34712(75) 


0.260 


5.35000(81) 


0.270 


5.35255(91) 


0.280 


5.35510(59) 


0.290 


5.35710(70) 


0.300 


5.35970(21) 


0.310 


5.36307(62) 


0.320 


5.36622(37) 


0.327 


5.36956(63) 


1/3 


5.37067(75) 



At first, we have considered interpolations with polyno- 
mials up to order /j,^ (see Table IIVI 1st to 4th row, for 
a summary of the resulting fit parameters) . We can see 
that data at /i^ < are precise enough to be sensitive to 
terms beyond the order /i^; indeed, a good x^/d.o.f. is 
not achieved before including terms up to the order /i^. 

As in the case of isospin chemical potential, we have 
performed a "constrained" fit. The largest interval for 
which a linear fit in /i^ works well is [(/i/ (TrT))^;^^, 0] with 
(/i/ (7rT))^;jj = —0.310^; we notice that such interval was 
larger {{fi / {nT))'^^^^ — —0.375^) in the case of an isospin 
chemical potential. Then, all available data were fitted 
by a 6th-order polynomial, with the constant term and 
the quadratic coefficient fixed at 5.32283 and —0.410, re- 



spectively (see Table IIVI 5th row) . 

Then, we have considered interpolations with ratios 
of polynomials: the one with the least number of pa- 
rameters for which we got a good fit is the ratio of a 
4th- to 2nd-ordcr polynomial (see Table HVl 6th row, and 
Fig.Hleft)). 

Finally, we have also tried the "physical" fit, as in the 
previous subsection, obtaining the following results for 
the "physical" fit ratio, Eq. ([3]): 

/3c(0) = 5.32373(90) , A = 8.140(32) , 

5 = 6.59(26), (7 = 7.201(35), (6) 

with x^/d.o.f.=0.51, see Fig. [Upright ) for a comparison of 
the fit to data for /3c (/i^)- 

Such interpolation permits, in principle, an extrapola- 
tion to real chemical potentials down to T = 0; from the 
parameters given in ^ one can get the extrapolation 
at T = of the pseudocritical quark chemical poten- 
tial: = TT^/C/B = 3.284(65) Tc(0). This result agrees 
within errors with the analogous one obtained in Ref. [l^ 
for SU(3) with n/ = 2 Wilson fermions on a smaller lat- 
tice and with smaller statistics, which turned out to be 
2.73(58) Tc(0). 

However, a comparison of different extrapolations 
shows that systematic effects become important well be- 
fore one approaches the T = axis. In Fig. [5] we have 
plotted the extrapolations to the interval < ^/{ttT) < 
0.5 of the following fits: 

• 6th-order constrained polynomial (5th row in Ta- 
ble [IV|; 

• ratio (4,2) of polynomials (last row in Table ITVl) : 

• "physical" fit ratio, Eqs. dH). 

The three curves agree as long as ^/{ttT) < 0.1, but 
then the 6th-order constrained polynomial deviates from 
the other two curves. Therefore results extrapolated to 
larger values of /i are not reliable. One could take the 
analogous results obtained at finite isospin chemical po- 
tential as a guiding reference, concluding that Fade like 
fits are to be preferred; however one cannot exclude that 
such argument may be wrong because of possible system- 
atic differences between QCD at finite quark and isospin 
chemical potentials. 

III. COMPARISON OF THE CURVATURES OF 
THE CRITICAL LINES 

In the present Section we will focus on the curvature 
of the pseudocritical line at /i = 0, which is the quan- 
tity with the least ambiguity related to the procedure 
of analytic continuation and for which a clear agreement 
among the determinations obtained by various different 
methods has been shown in previous literature [l^, . 
Our purpose is to determine how it changes when switch- 
ing from a theory at finite quark chemical potential to a 
theory at finite isospin chemical potential. 
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TABLE IV: Parameters of the fits to tfie pseudocritical couplings in finite density SU(3) witii Uf = 2 on a 16^ x 4 lattice 
with fermionic mass am= 0.05, according to the fit function /3c(At^) = (ao + ai(^/(7rr))^ + a2(fJ-/{TvT)f + aai/J./ {TvT)f)/{l + 
a4{fi/(TTT))^). Blank columns stand for terms not included in the fit. The asterisk denotes a constrained parameter. Fits are 
performed in the interval (fJ-/ {t^T))^-^^,0]; the last column gives the value of (/i/(7rT))^ 



/ mm ■ 



1 

min 



ao 



ai 



a2 



xVd.o.f. 



-0.310^ 

-1/32 
-1/32 
-1/32 
-1/32 



5.32189(78) 

5.32283(83) 

5.3242(13) 

5.3226(12) 

5.32283* 

5.32394(98) 



-0.4262(90) 
-0.410(10) 
-0.314(44) 
-0.446(86) 
-0.410* 
25.736(24) 



0.92(35) 
-1.7(1.7) 
-0.76(13) 
-1.05(61) 



-14.4(9.7) 
-8.7(5.4) 



4.9002(45) 



2.87 
0.63 
0.85 
1.41 

0.65 
0.60 



5.38 



5.36 



P. 5.34 



5.32 



5.3 



fit ratio (4,2) 



lattice 16 x4 



Quark ciiemical 
potential 




0.1 



5.38 



5.36 



P. 5.34 



"physical" fit 
ratio 

3 

lattice 16 x4 



Quark chemical 
potential 



(|X/(7lT))^ 




5.32 



(^l/(7lT))^ 



FIG. 4: Fits to the pseudocritical couplings in finite density SU(3) with n/ = 2 on a 16"^ x 4 lattice with fermionic mass 
am=0.05: ratio of a 4th- to 2nd-order polynomial (left) and "physical" fit according to the function Q (right). The dashed 
vertical line indicates the boundary of the first RW sector, Im(/i)/(7rr) = 1/3. 
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FIG. 5: Extrapolation to real quark chemical potentials of the 
6th-order constrained, ratio (4,2) of polynomials and "phys- 
ical" ratio fits (only the borders of the 95% CL band have 
been reported). 



In particular we want to determine the dependence of 
the critical temperature Tc{fiq, fiiso) at the quadratic or- 
der in /iq and /Xiso, which is determined by the two cur- 
vatures alone. Indeed one can show that, for two degen- 
erate flavors, the theory must be even under reflection of 
liq and /Ziso separately (see Ref. [ll| for a discussion on 
this point), so that the mixed /Zq/ijso term is absent and 

Tciflq, fiiso) = rc(0)+^q/4+Aso/4o + C(Mg/iso'AigAiL) ■ 

(7) 

In order to determine the two curvatures and compare 
them in a consistent way, we have performed a common 
fit to all the pseudocritical couplings at imaginary po- 
tentials reported in Tables U and IIIII with the following 
function 

/3c(m,, f^iso) = /3c(0) + aq (g;)' -f ai,o (^) ' (8) 

including as many data points, both at imaginary and 
real (when available) chemical potentials, as compatible 
with a reasonable value of x^/d.o.f. As a matter of fact 
the ranges of included chemical potentials coincide with 
those for which a linear fit works well separately for /i^ 
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or (see the second row of Tables |TT] and I IV I) , the 
only difference in this case being that /3c (0) is taken as a 
common parameter. The results of the fit are 



aq = -0.3997(87) 
^c(O) = 5.32370(57) , 



fli.o = -0.3606(67) 
xVd.o.f. = 0.93 . (9) 



We notice that a, and a 



are not compatible within 
errors and deviate from each other by about Aa. 

A convenient way to report the two curvatures is in 
term of dimensionless quantities, as follows: 



-^c(M(/5 Miso) 



1 



The parameters Rq and i?iso can be obtained respectively 
from aq and aiso, in particular one has 



'g/iso 



1 d a 



c(0) 



^q/iso 



Nr 



2/3c(0)3/3i(/3c(0),TO, 



(11) 



where a is the lattice spacing and Pl — o-idgo/da) is 
the lattice beta-function. Making use of the perturbative 
two- loop expression for we get 



Rq = -0.515(11) , i?iso = -0.465(9) 



(12) 



It is interesting to compare our results with those of 
previous studies. In Ref. 44] the same discretization 
and bare quark mass have been adopted for QCD at 
real isospin chemical potential; their result, when re- 
ported in the same units as ours, is i?iso = 0.426(19): the 
marginal discrepancy can be explained in terms of either 
the inexact R-algorithm or the smaller spatial volume 
used in Ref. Q. Rq = -0.500(34) has been obtained in 
Ref. [5)] for the same theory with a smaller fermion mass, 
am = 0.025: this is compatible with our result, show- 
ing that Rq has mild dependence on the quark mass. In 
Ref. 14] a value Rq = —0.38(12) has been reported mak- 
ing use oi rif = 2 Wilson fermions: the agreement, even 
if within quite large errors, is encouraging if we consider 
the completely different fermion discretization. Instead, 
as it is well known, the curvature changes significantly if 
we change the number of flavors; for instance for Uf = 4 
QCD one obtains Rq = -0.792(10) 

Our determinations of Rq and i?iso are clearly affected 
by the systematic error related to the choice of the two- 
loop expression for anyway such error disappears if 
we consider the ratio 



R. 



g — ISO 



Rq — Rh 

Rn 



= 0.098(26) , (13) 



which we consider as our final estimate for the difference 
in the curvature of the critical line between the theory 
at finite baryon density and the theory at finite isospin 
density. In order to appreciate the difference, in Fig. [5] 



we report the corresponding linear extrapolations to real 
chemical potentials. 

In previous studies the two curvatures revealed to 
be compatible within errors [13, El]. This is also the 
expectation in the limit of a large number of colors 
Nc 0, 113, 113 ■ indeed the two curvatures are ex- 
pected to be the same at the leading order 1/Nc [11] 
(the curvature itself is expected to vanish as Nc — >■ cxd). 
Therefore, we can consider the deviation that we find as 
the first evidence for an 0{1/N^) difference between the 
two theories at small chemical potentials. Rq-iso-, being 
the ratio of an 0{1/N^) to an 0{1/Nc) quantity, is ex- 
pected to be 0{1/Nc): this is compatible with the fact 
that it turns out to be of the order of 10%. It would 
be interesting to explore how results change for different 
values of Nr- 



5.32 




5.28 



5.26 



5.24 



5.22 - 



— isospin 
-- quark 
• simulations at real isospin 



J I I L. 



0.1 0.2 0.3 0.4 

^1/(7IT) , ^l. /(TlT) 



0.5 



FIG. 6: Comparison between the extrapolations to real quark 
and isospin chemical potential of the fits linear in ^/(tvT)'^. 
Data points (circles) are the results of Monte Carlo simula- 
tions performed directly at real isospin chemical potential. 



IV. ORDER OF THE PHASE TRANSITION 

As already stressed in the Introduction, the nature of 
the pseudocritical line at imaginary /ig, may be strongly 
influenced by the order of the RW endpoint [H, [l^, [H, 
[29| . i.e. the point at which the RW line taking place in 
the high-T region for Im(/ig)/T = 7r/3 meets the analytic 
continuation of the physical pseudocritical line. If the 
endpoint is first order then it is actually a triple point 
and at least the part of the pseudocritical line which is 
closest to the endpoint is expected to be first order. 

In the case oi Uf — 2, with the same regularization 
and temporal size [Nt — 4) used in the present study, it 
is known that the RW endpoint is first order, in the low 
mass region, for am < amn with amn = 0.043(5) f28j . 
That means that the mass used in the present work. 
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"L= 12 



i>-on/(jiT) = 0.425 
»-> n/(JiT) = 0.300 
&-nH/(jiT) = 0.165 i 




FIG. 7: Normalized plaquette distributions at the pseudo- 
critical coupling for diflterent values of the isospin chemical 
potential. 




FIG. 8: Normalized plaquette distributions at the pseudocrit- 
ical coupling for different spatial lattice sizes and ^iso/{T^T) = 
0.475i. 



am = 0.05, is close to the tricritical value but slightly 
on the second order side, so we do not expect the an- 
alytic continuation of the pseudocritical line to become 
first order as we approach the RW endpoint. This is com- 
patible with the fact that we have not observed signals of 
metastable behavior or double peak distributions along 
the line; only a strengthening of the transition can be seen 
as the RW endpoint is approached, as a consequence of 
the closeness of the tricritical point. 

If one conjectures that an imaginary fiiso may 
strengthen the transition in the same way as an imagi- 
nary does, then, since the range available for Im(/iiso) 
is larger than that available for Im(/iq), one may ex- 
pect that a first order transition could be manifest at 
some stage along the pseudocritical line at imaginary 
^iso- Such conjecture is well founded, since simulations 
at real isospin chemical potential have shown that indeed 
the effect of small positive values of /i?^ is a weakening 
of the transition 50, 51]. 

In order to explore this possibility, we have reported 
in Fig. [7] the plaquette distributions at the pseudocritical 
coupling for a few different values (both real and imag- 
inary) of ^iso/(7rT). It is evident that for the largest 
values of /iiso a double peak structure develops, hinting 
at the presence of a first order transition. 

In order to confirm that by a finite size scaling analy- 
sis, we have repeated simulations for the largest value of 
Im(/Xiso), /iiso/(7i'T) = 0.475i, on two other lattice sizes, 
L = 12 and L — 20. Both the scaling of distributions 
and the scaling of susceptibilities confirm the first order 
nature of the transition for this value of ^iso : the well in 
the double peak distribution of the plaquette deepens as 
L increases as expected (see Fig. and the maxima of 
the plaquette susceptibility scale linearly with the spatial 
volume (see Fig. IH]). 

Therefore we conclude that, for the present discretiza- 
tion and value of the quark mass, the transition is surely 
first order at /iiso/(7rT) = 0.475i and there is possibly 



a critical point along the line at some smaller value of 
Im(^iso). Such non-trivial behavior resembles what hap- 
pens for quark chemical potentials [1^ 113, [28|, (2^ and 
may have consequences on the general structure of the 
QCD phase diagram which should be further investigated 
in the future. 



0.4- - 

0.3- 

0.2- 

0.1- 

qI , \ , \ , \ , \ 

2000 4000 6000 8000 

FIG. 9: Maxima of the plaquette susceptibility as a function 
of the spatial volume for /iiso/ (tT) = 0.475i. 

V. CONCLUSIONS 

In this paper we have considered QCD with two de- 
generate fiavors of bare mass am — 0.05, corresponding 
to a pion mass TOtt ^ 400 MeV, on a 16"^ x 4 lattice, at 
nonzero quark or isospin density. 

Our investigation developed along three main lines: 

• localization of the pseudocritical line in the 
temperature-chemical potential plane, for the two 
cases of quark and isospin density; 
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• comparison of the curvatures of the two critical 
hues at the point of zero chemical potential; 

• study of the order of the phase transition along the 
two critical lines. 

To study the localization of the critical line we adopted 
the method of analytic continuation. The result we got 
is similar to what we found in the case of four degenerate 
flavors and is common to the case of quark and isospin 
density: deviations from the linear behavior in fj? of the 
critical lines are clearly seen for /x^ < and are nicely 
described by several analytic functions. However, the ex- 
trapolations to positive /i^ overlap, within errors, only as 
long as Aiiso/(7rT) ~ 0.2 and Hq/inT) ~ 0.1. The com- 
parison with direct numerical simulations performed at 
real isospin chemical potentials leads to a preference for 
extrapolations based on Pade approximants; such sug- 
gestion, which is agreement with our previous studies, 
could be taken as a guiding principle also in the case of 
nonzero quark density. 

We have performed a careful determination of the cur- 
vatures of the two critical lines at zero chemical potential. 
In order to compare them in a consistent way, we have 
performed a common fit of the pseudocritical couplings, 
taking the critical /? at zero chemical potential as a com- 
mon parameter. We have found that the curvature of the 
isospin critical line is larger than that of the quark criti- 
cal line by about 4cr, the relative difference being about 
10%. The outcome of previous studies |3, Ell 
favour of a substantial agreement between the two cur- 
vatures, as expected in the limit of a large number of 
colors Nc l4 6| . [ 47l. IZsl . l49j and in particular at the leading 
1/Nc order [ia]. The deviation that we find is therefore a 
first evidence for an 0{1/N^) difference between the two 
theories at small chemical potentials. The order of mag- 
nitude of the relative deviation Rq^iso, which is the ratio 
of an 0{l/N^) to an 0{l/Nc) quantity (the curvature 
itself) is compatible with it being an 0{1/Nc) quantity. 
It would be interesting to explore how results change for 



different values of Nc- 

Finally, we have studied the order of the transitions 
along the critical lines. For the case of nonzero quark 
density, we have found no clear signatures of a first order 
transition, in agreement with the expectations after re- 
cent findings in the literature and in consideration of the 
quark mass adopted in this work, which is larger than 
the tricritical mass found in Ref. t28j. The only effect 
we could see was a strengthening of the transition when 
the RW point is approached at imaginary quark chemical 
potential. 

A phenomenon emerged from our investigation is in- 
stead that an imaginary isospin chemical potential can 
strengthen the transition, similarly to what happens for 
a quark chemical potential and in agreement with the 
fact that simulations at real isospin chemical isospin have 
shown that small positive values of /i? ^ weaken the tran- 
sition [sO, Hlj. Moreover we have found clear evidence 
that, in this case, the transition becomes first order for 
large enough imaginary chemical potentials, but before 
one reaches the RW-like transition which is found, in the 
high-T region, for ^iso/{TrT) ~ 0.5 ^l|, implying a possi- 
ble second order critical point along the line. Such behav- 
ior is very similar to what found for a quark chemical po- 
tential, hinting at a common underlying physical mech- 
anism, and could have various non-trivial consequences, 
on the shape of the critical line and on the general struc- 
ture of the QCD phase diagram, that should be further 
investigated in future studies. 
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